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We prove that the inclusive single-gluon production cross section for a hadron colUding with a 
high-density target factorizes into the gluon distribution function of the projectile, defined as usual 
within the DGLAP coUinear approximation, times the cross section for scattering of a single gluon 
on the strong classical color field of the target. We then derive the gluon-proton (nucleus) inelastic 
cross section and show that it is (up to logarithms) infrared safe and that it grows slowly with center 
of mass energy. Furthermore, we discuss jet transverse momentum broadening for the case of nuclear 
targets. We show that in the saturation regime, in contrast to the perturbative regime, the width of 
. . . the transverse momentum distribution is infrared finite and grows rapidly with energy and rapidity. 

CN ! In both regimes, however, transverse momentum broadening exibits the same A dependence. 

O 

o 

^ , i. introduction 

Oh; 

' Understanding the behavior of hadronic cross sections at very high energy is one of the major unresolved problems 
in QCD. Even though Regge theory can, in principle, predict the energy dependence of hadronic cross sections and 
there are many phenomenological models, based on Regge theory, which are somewhat successful in describing the 
■ data, the relation of Regge theory to QCD is not well understood. Therefore, it would be very important to be able 
, to calculate the high energy behavior of hadronic cross section from QCD itself. However, it is believed that total 
^ ■ cross sections are intrinsically non-perturbative and not amenable to perturbative QCD methods [Q. In this note, 
— , we consider the simpler problem of (real) gluon-proton total inelastic cross section in the very high energy (small x) 
' limit, using the effective action and classical field method developed recently, and show that its growth with energy 
CO I is inhibited as compared with that expected from perturbative QCD. 

At very small x, a hadron is a Color Glass Condensate due to the condensation of gluons into a coherent state with 
characteristic momentum of Qs{x) In other words, most of the gluons in the wave function of a hadron have 

momenta of order Qs{x). As we go to higher energies, Qs{x) increases and eventually will become much larger than 
^QCD so that as{Qs) <C 1 and weak coupling methods can be used. Even though the theory may be weakly coupled, 
_ . it is still non-perturbative in the sense that one has high gluon densities and strong classical color fields associated with 
II ■ them so that the standard perturbative QCD breaks down. Much progress has been made in developing a formalism 
CIh' which describes this weakly coupled, though non-perturbative region of QCD. It generalizes the standard coUinear 
^ factorized (leading twist) expressions for hadronic cross sections and allows one to calculate hadronic cross sections 
in an environment where higher twist (high gluon density at small x) effects are important. We refer the reader to 
and references therein for a review of this formalism. 

One can use this classical field method to calculate gluon production at high energy Q. In this Letter, we use 
the results of ^ to prove (coUinear) factorization of the inclusive cross-section for a "dilute" hadron impinging on a 
' dense target. Based on that result we then derive an expression for the gluon-proton (nucleus) inelastic cross section 
and for its energy dependence at very high energies. Using our results for the gluon-proton (nucleus) inelastic cross 
section, we consider the transverse momentum broadening of gluons due to scattering from the strong classical field 
of a nucleus and show that it is infared finite, energy dependent and scales like A^^^. 

II. GLUON-PROTON INCLUSIVE CROSS SECTION 

In the classical field and effective action approach to hadronic (nuclear) collisions at high energy, one solves the 
classical Yang-Mills equations of motion in the presence of random color charges created by the valence quarks and 
gluons at high x. One then averages over these color charges with a Gaussian weight to compute physical quantities. 
The classical fields of the colliding hadrons (nuclei) before the collision are given by the single-hadron (nucleus) 
solutions which, in light cone gauge, are 

^1,2 = -f/i,2(xi)5'f//,2(a^±) • (1) 
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These fields serve as the initial conditions for solving the Yang-Mills equations of motion in the forward light cone; 
they were solved in Q for the case of asymmetric collisions where the classical field of one of the colliding sources is 
much stronger than the classical field of the other source. 




proton (Q^> p^) 



FIG. 1. Gluon production in inclusive pp scattering at rapidity far from the target proton (i.e., the strong color field). 
The produced gluon field in the forward light cone region is given by 



A^{t,x^) = ±x^U{xi_)l3{T,Xi_)U\x^) . 



(2) 



We have chosen the gauge condition x+A^ + .t^A+ = 0. Here, r — \/2x+x^ denotes proper time and x±^ is the 
transverse coordinate. The ?7's are rotation matrices in color space, to be specified shortly. At asymptotic times, 
T ^ oo, the fields (3 and /3* are given by superpositions of plane wave solutions, 



(3(t CO, x±) = 
j3'{T oo, x±) = 



1 



(2^)' V2[J^ 



(27r)2 uj 



{ai{p^)e'P^-^^-'^^ + C.C.} , 
{a2ip±)e'P^-''^-"^^ + c.c.} 



The number distribution of produced gluons at rapidity y and transverse momentum p± is given by 



dN 



(Pp±dy (27r) 



rtr 



(|ai(p±)|V|a2(p±)f 



(3) 
(4) 

(5) 



One can obtain an analytical solution of the classical Yang-Mills equations if one of the fields (the "projectile") is 
much weaker than the other (the "target"). This situation is realizedphysically when \y — yt\ 3> \y ~yp\i because of 
the renormalization group evolution of the gluon density in rapidity l5|J^,y|. In that case, it turns out ||^ that the J7's 
appearing in eq. (^ are just the C/2's from eq. (|l]); that is, to leading order in the weak field, the plane waves in the 
forward light cone are just gauge rotated by the strong field 



U2ix±,y) ^T'exp -ig / dy'<i>2ixA_,y') 



(6) 



Vt 



Here, we assumed that the target is moving along the negative z-axis, i.e. at negative rapidity yt < 0. In the Color 
Glass Condensate model, we now have to average the squared amplitudes over the gauge potentials tf>i,2(a;j_, y) using 
a Gaussian weight 



|ai,2|')= y"l?$il?$2|ai,2($i,$2)|' 

tr (Vi<fi(xi,y'))' 



X exp 



- dy' d^xi_ 



dy' / d'^x± 



tr {S/l^2{x ±,y')y 
9'^fJ'l{x_L,y') 



(7) 
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Here, fj, denotes the density of color charge in the sources per unit of transverse area and rapidity. The radiation 
number distribution (|^) turns out to depend only on the integrated color charge densities of the sources, 

Vp y 
Xi{y,pl) = I dy'fil{y',pl) , X2{y,pl) = J dy'fil{y',pl) . (8) 

y yt 

Due to the evolution in rapidity, Xiiv) X^iu) if \y — yt\ ^ |y — Vpl- In terms of the gluon distribution function in 
the projectile or target proton respectively. 

Nr. 1 } , .2 



xiAy:P±) = jpTZTi:^ J '^^9pAx,p±) , (9) 

Xo 

where Xo is on the order of p± cosh(y)/-ys. In the weak-projectile limit the averaging (|^) over the gauge potentials 
with Gaussian weight can be done Q . At high transverse momentum, p± ^ Qt (wh ere Qt is the saturation momentum 
of the target), one recovers the standard result from perturbation theory l7| Jl0| , [ll|] , 

- s c s c I ± I dx'gp{x\kl) dx"gt{x",{p^-k^f) . (10) 



d'^p±dy p\ J 7r2 k^ip^-ki:, 

This is to be expected, of course, because the gluon occupation numbers in the target proton become small at p±^ ^ Qt, 
i.e. that field becomes weak as well and so can be treated perturbatively. 

In the coUinear limit {kj_/p±^ 0), the first factor in ( p^ is the DGLAP splitting function for gluons. It 
evolves the gluon distribution function of the projectile from the scale fc^ to the scale p\, 

p± 

~;{2~^ J dxgpix^kl) . (11) 



Thus, in the coUincar limit cq. ( pl]| ) simply turns into 

dN N^ 1 



xgp{x,pl) / dx"gt{x",pl) . (12) 



d'^pjLdy 7ri?2 - 1 p4 

Integrating over the impact parameter space (simply a factor of ttR^) and dividing by the flux of incoming gluons in 
the projectile proton at the hard scale p^, which is simply xgp{x,p'^), we obtain the inclusive gluon-proton differential 
cross section at large momentum transfer 

Xo 

This is the standard expression for the gluon-proton inclusive cross section [^o|jlJ] in the one-gluon exchange ap- 
proximation which diverges like l/i^ = l/p\ at small momentum transfer, cr^p'* grows like a power of energy if 
xg{x) ^ l/x^ with S > 0. 

Next, we consider the radiation distribution (|^) at small transverse momentum, < Qt, but large enough so that 
the projectile proton is still in the weak-field regime p± '3> Aqcd- In this regime m 



dN 1 1 d^k± p2 ^ 



^asNc / — —JTI i-T2 / dxgp{x,k^) (14) 



d^p±dy Att p\ J TT^ k\{p^--ki_Y 

xo 
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- ^:^^9pix,P±) ■ (15) 

This result is valid to leading order in a^xi, but to all orders in 0^X2- The second line applies in the limit of nearly 
collinear splitting {k±/pj_ — > 0), where we can employ eq. (|ll|). Thus, the inclusive gluon production again factorizes 
into the gluon distribution function of the projectile at the scale times the cross section for scattering of a gluon 
on the high-density target. 
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To obtain the gluon-proton cross section, we again integrate over the impact parameter and divide by the flux of 
incoming coUinear gluons in the projectile at the scale p^, which is xgp{x,p\). We get 

— ^ = 1 T^R 16 



In the saturation regime, the cross section is of order 1 rather than as in (13) because the occupation number of the 
target ^ l/ag cancels one power of the coupling in both the amplitude and the complex conjugate amplitude. Note 
that the above result holds for p± ^ Age d only. To get the gluon-proton total inclusive cross section, we integrate 
(H) which gives 



alf ^ -nR' log (QV^Icd) + 0{al) . (17) 



Parameterization of the saturation momentum like a power of energy, ~ s'^ as done in |T^] then leads to a 
logarithmically growing cross section 

fjgp ttR^ logs . (18) 



On the other hand, if we use a DLA DGLAP type parameterization Qt ~ exp \J\og l/x, the cross section would grow 
like a square root of energy. Therefore, it is clear that different parameterization of the target saturation scale Qt 
would lead to a different growth of the cross section with energy. Curiously enough, assuming Qt to grow like 5^°^ ^ 
would lead to growth of the gluon-proton cross section like log^ s. In order to determine the energy dependence of 
Qt rigorously, one would need to formally define it in terms of a gluonic two-point function and solve the non-linear 
evolution equation that Qt follows. This is however beyond the scope of this work. 

Strictly speaking, our results are correct for the cross section per unit area, da/dPb. What we have shown here is that 
the growth of this cross section at a fixed impact parameter is inhibited due to high gluon density effects. In principle, 
as we go to larger impact parameters 6, the gluon density becomes smaller and smaller until our classical formalism 
reduces to the standard pQCD where the gluon density is not that large. Since large impact parameters (where the 
gluon density is low and so the classical approach is not valid) are believed to give the dominant contribution to 
the total cross section, one should understand the area appearing in ([l^ ) as a parameter which we can not precisely 
determine. A genuine non-perturbative (strong coupling vs. our weak coupling methods) calculation would presumably 
enable one to determine this factor and its energy dependence. 

The above equations hold equally well for a nuclear target. In eq. ( p^ ) one just has to replace the gluon distribution 
function of the target proton by that for the nucleus; in the absence of shadowing ^,|l^, gA[x, Q^) — Ag[x, Q^)- On 
the other hand, in eqs. ([l^ , p^ one substitutes the radius of the proton by that of the nucleus, Ra — A^^^R. 

It is now straightforward to compute the transverse momentum broadening of the incoming gluon jet traversing 
the target nucleus |l4yi7f| . The transverse momentum broadening is given by fl^ 

= (19) 
where (tA{b)] is the nuclear thickness function averaged over impact parameters: 



tA{h)) = ^=pL, (20) 



A 

3 



with L the average thickness of the nucleus, and with p ~ 0.15 fm the density of nucleons in the nucleus. 
In perturbation theory, using ( p^ ) and DGLAP evolution. 



asNc f , , 9 ^ <i_ 

rflogpi 



dxgA{x,p\) = ^i„„^2 ^9a{x,p^^ , (21) 



one obtains 

{p\) = 47r2 a, -^^i^ {xg{x, Ql,a,) - xgix, Ql,,J] pL . (22) 



Thus, (p^) grows proportional to A^^^ if nuclear shadowing is disregarded p^ , p7t . Appearance of xg{x, Q^i„) signifies 
sensitivity to physics at small momentum transfer. For a discussion of what value of x is to be used in (^), see []l^ , p7t . 
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Transverse momentum broadening in the saturation regime is quite different. Using (jl^) and ( p^|j2C| ), we find that 

(pl)o,9l+o{a.), (23) 

where Q\ denotes the nuclear saturation scale which, from the definition of X2 ~ Qa/^^s ®' A-"^/^ times larger 
than that for a proton. Thus, in the saturation regime jet broadening grows with the same power of the nuclear mass 
number A as in the perturbative regime ~ A^^^. Also, our derivation shows that the jet transverse momentum 
broadening is energy dependent even though we can not determine its energy dependence in a model independent 
way. However, for Q\ ~ s'^ p5| , we obtain a strong power-law increase of transverse momentum broadening. We 
would like to emphasize jl^ that since the saturation scale of the target is expected to be much larger in the 
forward rapidity (projectile fragmentation) region, the jet transverse momentum broadening will be much larger in 
the forward region than in the central rapidity region, contrary to the prediction (|2^ ) from perturbation theory. 



III. CONCLUSION 



In summary, we have shown that inclusive gluon production from a hadron scattering off a high-density target 
(with "saturated" gluon density) factorizes into the gluon distribution of the projectile times the cross section of the 
beam of incoming coUinear gluons on the dense target. We derived the gluon-proton inclusive cross section in the high 
energy limit. We have shown that the cross section grows, with reasonable, HERA compatible parameterizations 
of the saturation momentum, only logarithmically with energy rather than power like as expected from perturbation 
theory. This "unitarization" of the cross section is due to the strong classical fields of the target generated by the high 
gluon density (higher twist) effects. We have also considered the transverse momentum broadening of the gluon jet 
passing through a nuclear target. We have shown that it scales like A^^^ in both perturbative and saturation regimes 
and that it is infrared finite. We predict that the jet transverse momentum broadening will be larger in the forward 
rapidity region and that it increases with energy. 
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